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ABSTRACT 

High-resolution optical spectra of eighteen blue horizontal-branch (BHB) stars in 
the globular cluster M15 indicate that their stellar rotation rates and photospheric 
compositions vary strongly as a function of effective temperature. Among the cooler 
stars in the sample, at T e g- ~ 8500 K, metal abundances are in rough agreement with 
the canonical cluster metallicity, and the vsini rotations appear to have a bimodal 
distribution, with eight stars at vsini < lSkms" 1 and two stars at vsini ~ 35kms _1 . 
Most of the stars at T e g > 10000 K, however, are slowly rotating, vsini < 7kms~ 1 , 
and their iron and titanium are enhanced by a factor of 300 to solar abundance levels. 
Magnesium maintains a nearly constant abundance over the entire range of T e g, and 
helium is depleted by factors of 10 to 30 in three of the hotter stars. Diffusion effects 
in the stellar atmospheres are the most likely explanation for these large differences in 
composition. Our results are qualitatively very similar to those previously reported 
for M13 and NGC 6752, but with even larger enhancement amplitudes, presumably 
due to the increased efficiency of radiative levitation at lower intrinsic [Fe /H] . We also 
see evidence for faster stellar rotation explicitly preventing the onset of the diffusion 
mechanisms among a subset of the hotter stars. 

Subject headings: globular clusters: general, globular clusters: individual (NGC 7078), 
stars: horizontal-branch, stars: abundances, stars: rotation 
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1. Introduction 

M15 (NGC 7078) is one of the most metal-poor globulars known, with a metallicity 
[Fe/H] ~ —2.4 dex measured from red giant abundances (Cohen 1979, Sneden et al. 1997). Like 
many other such clusters, M15's horizontal branch lies predominantly bluewards of the instability 
strip, and color-magnitude diagrams (Buonanno et al. 1983; Durrell & Harris 1993) show an 
extended "blue tail" reaching T e g as high as 20000 K, which is separated from the horizontal part 
of the HB by a "gap" in the distribution of stars along the HB (Moehler et al. 1995). Similar gaps 
appear in the CMDs of M13, M80, NGC 6752, NGC 288, and other clusters, but are difficult to 
explain via standard models of RGB mass loss or HB evolution. 

Recently, attention has focused on atmospheric effects as a possible explanation for these 
photometric features (Caloi 1999; Grundahl et al. 1998). At T e g ~ 10000 K, they suggest, the 
stellar atmospheres become susceptible to diffusion effects, and thus develop chemical peculiarities 
similar to those that appear in main-sequence CP stars. The resulting changes in atmospheric 
opacity alter the emitted flux distributions and thus the measured photometry of the hotter stars, 
giving rise to the gaps. These claims have been bolstered by measurements of large photospheric 
abundance anomalies among hotter BHB stars in M13 (Behr et al. 1999) and NGC 6752 (Moehler 
et al. 1999), which show 30 to 50 times the iron abundance expected for these metal-poor clusters. 
In M13, the transition from normal-metallicity cooler stars to metal-enhanced hotter stars is 
remarkably abrupt, and coincides with the location of the gap, further supporting the surface-effect 
hypothesis. 

It appears, however, that stellar rotation also plays some role in this process. Theoretical 
treatments of the diffusion mechanisms (Michaud et al. 1983) suggest that circulation currents 
induced by higher rotation velocities can easily prevent abundance variations from appearing. 
This prediction is borne out by measurements of vsini for the M13 stars (Behr et al. 2000), which 
show that although the cooler stars exhibit a range of vsini, some as high as 40kms _1 , the hot 
metal-enhanced stars all show very low rotations, vsini < lOkms -1 . This correlation suggests 
that slow rotation may be required in order for the metal enhancement and helium depletion 
appear in the photosphere. 

Although these results from M13 and NGC 6752 imply that we are on the right track 
towards explaining the photometric peculiarities, BHB stars in many other clusters will have to be 
analyzed in a similar fashion before we can make any firm claims. In particular, since the radiative 
levitation that causes the observed metal enhancements is thought to depend strongly on the 
intrinsic metallicity of the atmosphere, we should study clusters spanning a range of [Fe/H], to see 
whether the onset and magnitude of the enhancements vary. To this end, we have observed BHB 
stars in five clusters in addition to M13. In this Letter, we describe the results for both rotation 
and photospheric abundances in M15. The data for the other clusters — M92, M68, NGC 288, 
and M3 — are being analyzed currently, and will be presented in a future publication. 
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2. Observations and Reduction 

The eighteen stars in our sample were selected from Buonanno et al. (1983), and are listed in 
Table 1. We acquired supplementary Stromgren photometry of the target stars using the Palomar 
60-inch telescope, in order to better constrain the effective temeratures. The stars generally lie in 
the cluster outskirts, where crowding and confusion are less of a problem than towards the core, 
and our CCD imaging confirmed an absence of faint companions. 

The spectra were collected using the HIRES spectrograph (Vogt et al. 1994) on the Keck I 
telescope, during four observing runs on 1997 August 01-03, 1997 August 26-27, 1998 June 27, 
and 1999 August 14-17. A 0.86-arcsec slit width yielded R = 45000 (v = 6.7 km s" 1 ) per 3-pixel 
resolution element. We limited frame exposure times to 1200 seconds, to minimize susceptibility 
to cosmic ray accumulation, and then coadded four frames per star. S/N ratios were on the order 
of 30 — 60 per resolution element. 

We used a suite of routines developed by J.K. McCarthy (1988) for the FIGARO data 
analysis package (Shortridge 1988) to reduce the HIRES echellograms to 1-dimensional spectra. 
Frames were bias-subtracted, flat-fielded against exposures of HIRES' internal quartz incandescent 
lamps (thereby removing much of the blaze profile from each order), cleaned of cosmic ray hits, 
and coadded. A thorium-argon arc lamp provided wavelength calibration. Sky background was 
negligible, and 1-D spectra were extracted via simple pixel summation. A lOth-order polynomial 
fit to line-free continuum regions completed the normalization of the spectrum to unity. 

3. Analysis 

The resulting spectra show a few to over 140 metal absorption lines each, with the hottest 
stars showing the largest number of lines. Line broadening from stellar rotation is evident in a 
few of the stars, but even in the most extreme cases, the line profiles were close to Gaussian, so 
line equivalent widths (W\) were measured by least-square fitting of Gaussian profiles to the data. 
Equivalent widths as small as 10 mA were measured reliably, and errors in W\ (estimated from 
the fit y 2 ) were typically 5 mA or less. Lines were then matched to the atomic line lists of Kurucz 
& Bell (1995), and those that were attributed to a single species (i.e. unblended) were used to 
determine radial velocity v r for each of the stars. On the basis of v r , all of the targets appear to 
be cluster members. 

To derive photospheric parameters T e g and log g, we compared the published photometry and 
our own Stromgren indices to synthesized colors from ATLAS9 (Kurucz 1997), adopting a cluster 
reddening of E(B — V) = 0.09. Temperatures are well-determined (± a few hundred K) for the 
cooler stars, but are somewhat less firm (as much at ±1000 K) for the hotter stars. This situation 
will improve with a more sophisticated treatment of the various Stromgren colors, and by using 
transitions with different excitation potentials \ to constrain T c g. We estimated surface gravities 



- 5 - 



using the AB V flux method (Oke & Gunn 1984), which relates log g, stellar mass M#, distance d, 
and photospheric Eddington flux H u at 5480 A. Since ^45^(5480^4) = V magnitude, we can derive 

log g = 3.68 + log(M*/M ) + log(#„) — (M — m) v + 0.4Vb , 

where unextincted magnitude Vq = V + Ay = V + 0.30 for M15. We assumed = O.6M as 
a representative BHB star mass, used a distance modulus of 15.26 (Silbermann & Smith 1995), 
and drew the B. v values from the ATLAS9 model atmospheres (Kurucz 1997), iterating until log g 
converged. The resulting logg values agree well with model ZAHB tracks (Dorman 1993), except 
for the hotter stars, which are "overluminous" as described by Moehler et al. (1995). Table 1 lists 
the final photospheric parameters used for each of the target stars, as well as the heliocentric 
radial velocities. 

For the chemical abundance analyses, we use the LINFOR/LINFIT line formation analysis 
package, developed at Kiel, based on earlier codes by Baschek, Traving, and Holweger (1966), with 
subsequent modifications by M. Lemke. Our spectra of these very metal-poor stars are sufficiently 
uncrowded that we can simply compute abundances from equivalent widths, instead of performing 
a full spectral synthesis fit. Only lines attributed to a single chemical species were considered; 
potentially blended lines are ignored in this analysis. Upper bounds on [Ti/H] were determined for 
two of the hotter stars, which did not show any titanium lines, by assigning an equivalent width 
of 20 mA (double the strength of the weakest metal lines actually measured in those spectra) to 
the six strongest Ti II transitions and calculating the implied abundance. Microturbulent velocity 
£ was chosen such that the abundance derived for a single species (Fe II for most cases) was 
invariant with W\. For those stars with an insufficient number of lines to utilize this technique, we 
adopted a typical value of £ = 2kms _1 . We assumed a cluster metallicity of [Fe/H] = —2.4 dex 
in computing the initial model atmospheres, and then adjusted as necessary for those stars that 
turned out to be considerably more metal-rich, although these adjustments to the atmospheric 
input were found to have only modest effects (< 0.2 dex) on the abundances of individual elements. 

Measurement of v sin i broadening traditionally entails cross-correlation of the target spectrum 
with a rotational reference star of similar spectral type, but this approach assumes that the 
template star is truly at v sin i = 0, which is rare. Furthermore, given the abundance peculiarities 
that many of these BHB stars exhibit, it is difficult to find an appropriate spectral analog. Since 
we are able to resolve the line profiles of our stars, we instead chose to measure vsini by fitting 
the profiles directly, taking into account other non-rotational broadening mechanisms. We used 
bright unsaturated arc lines to construct an instrumental profile, and then combined it with an 
estimated thermal Doppler broadening of Skms" 1 FWHM and the previously-determined £. This 
profile was convolved with hemicircular rotation profiles for various vsini to create the final 
theoretical line profiles. Each observed line in a spectrum was fit to the theoretical profile using 
an iterative least-squares algorithm, solving for v sin i and line depth. The values for v sin i from 
different spectral lines generally agree quite well, once we removed helium lines and blended lines 
such as the Mg II 4481 triplet. 
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4. Results 

In Figure 1, abundance determinations for key chemical species are plotted as a function 
of stellar T e g. The values [X/H] represent logarithmic offsets from the solar values of Anders & 
Grevesse (1993). The error bars incorporate the scatter among multiple lines of the same species, 
plus the uncertainties in T e g, log g, £, W\ for each line, and [Fe/H] of the input atmosphere. 

In the cooler stars (T c g < 10000 K), the compositions are largely as expected. The iron 
abundance [Fe/H] averages —2.5, slightly below the value of —2.4 expected for this metal-poor 
cluster. Magnesium and titanium appear at [Mg/H] ~ —2.2 and [Ti/H] ~ —1.8, respectively, 
which are also reasonable for this environment. 

As we move to the hotter stars, however, the iron abundances change radically. Six of 
the eight stars at T e g- > 10000 K show solar iron abundances, [Fe/H] ~ 0, an enhancement of 
a factor of 300 from the cooler stars. (The other two hot stars show the same [Fe/H] as the 
cool stars, for reasons which will be discussed shortly.) Titanium, in a similar fashion, rises to 
[Ti/H] > 0, although there are hints of a monatonic increase with T e g- rather than an abrupt 
jump. Magnesium, on the other hand, appears to be unaltered, maintaining a roughly constant 
metal-poor level across the entire temperature range. The hotter stars also start to show helium 
lines, providing evidence of helium depletion, since [He/H] ~ at 11000 K but it then drops by 
factors of 10 to 30 for T cff = 12000-13000 K. 

Figure 1 also charts the values of v sin i derived for the target stars. Among the cooler stars, 
we find a range of rates, with most of the stars rotating at 15kms _1 or less, except for two stars 
at 29kms _1 and 36kms _1 , respectively. This sort of distribution of vsini is not what one would 
expect given a single intrinsic rotation speed v and random orientation of the rotation axes, since 
large sin % are more likely than small sini. Instead, the cool end of the HB appears to contain two 
rotational populations, one with v ~ 35kms _1 , and another with v ~ 15kms _1 , much like in M13 
(Peterson et al. 1995). 

For the hotter stars, there also appears to be a bimodal distribution in v, although the 
difference is less pronounced. Six of these eight stars show vsini < 7kms _1 , while the other two 
have vsini ~ 12kms _1 . Interestingly, these two faster-rotating hot stars are the same stars that 
show "normal" (metal-poor) iron abundances. 

5. Discussion 

These BHB stars in M15 exhibit abundance and rotation characteristics very similar to those 
reported for BHB stars in M13 (Behr et al. 1999 Sz 2000). The enhancements of metals except for 
magnesium, the depletion of helium, and the difference in maximum v sin i between hotter and 
cooler stars are shared by both clusters. The abundance anomalies in M15 are therefore likely to 
be due to the same diffusion processes that were invoked for the prior study — radiative levitation 
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of metals, and gravitational settling of helium, in the stable non-convective atmospheres of the 
hotter, higher-gravity stars, as hypothesized by Michaud et al. Unfortunately, the stars selected 
in M15 do not sample the immediate vicinity of the photometric gap as well as those in M13, so 
the association between the onset of diffusion-driven abundance variations and the gap is not as 
clear-cut, but the general trend still supports this association. 

There are, however, two notable differences between the results for M15 and those for M13. 
First, the magnitudes of the metal enhancements are somewhat different. Iron and titanium are 
each enhanced by ~ 2.5 dex in M15, while in M13, iron increases by only ~ 2dex, and titanium by 
~ 1.5 dex. Despite this difference, the enhancement mechanism yields the same final metallicities 
for the hot stars in both clusters: [Fe/H] ~ and [Ti/H] ~ 0. This correspondence suggests that 
the radiative levitation mechanism reaches equilibrium with gravity at or near solar metallicity, 
independent of the initial metal content of the atmosphere, as the metal lines become saturated 
and are thus unable to support further enhancements. 

Second, there is the issue of the two hotter M15 stars which do not show metal enhancement, 
denoted by circles in Figure 1. Their derived temperatures associate them with the hotter 
population, as does their photometry, which places them bluewards of the photometric gap. Their 
iron abudances, however, are < —2.5 dex, like the cooler unenhanced stars in the cluster, and 
stringent upper bounds on their titanium abundances again suggest that they are metal-poor. 
These two stars are also distinguished by having vsini ~ 12kms" 1 , nearly twice as large as any 
of the other hot stars. Although their temperatures and gravities are high enough to support 
radiative levitation, it appears that their faster rotations induce meridonal circulation, which keeps 
the atmosphere well-mixed and prevents the metal enhancements from appearing. Such sensitive 
dependence upon rotation speed was mentioned by Michaud et al, but these observations provide 
direct evidence that rotationally-driven mixing can and does directly influence the atmospheric 
composition. 

With these results from M15, we add to the growing body of evidence that element diffusion, 
regulated by stellar rotation, are at least partially responsible for the observed photometric 
morphology of globular cluster HBs. The findings from this cluster corroborate the prior work 
on M13 and NGC 6752, while adding some new twists which further illuminate the diffusion 
mechanisms. Analysis of many other clusters, spanning a range of metallicity and HB morphology, 
will be necessary before all of the implications of these effects can be known. 

These observations would not have been feasible without the HIRES spectrograph and the 
Keck I telescope. We are indebted to Jerry Nelson, Gerry Smith, Steve Vogt, and many others for 
making such marvelous machines, to the W. M. Keck Foundation for making it happen, and to a 
bevy of Keck observing assistants for making them work. J.G.C. and B.B.B. acknowledge support 
from NSF Grant AST-9819614. 
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Fig. 1. — Fe, Ti, Mg, and He abundances, and projected rotation velocity vsini for BHB stars over 
a range of T e fj. Abundances are plotted as log offsets from the solar abundances. Circles indicate 
the two hot fast-rotating non-metal-enhanced stars discussed in the text, and the open symbols in 
the [Ti/H] plot are upper bounds on the titianium abundance for these two stars. 
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Table 1. Parameters for program stars 



Star 


V 


B-V 


-^lincs 


V r 

(kms- 1 ) 


T e s 
(K) 


log9 


(kms- 1 ) 


v sin i 
(kms- 1 ) 


B124 


15.91 


0.15 


47 


-106.0 


7900 ± 100 


2.80 ±0.05 


2 


6.38 ± 0.25 


B558 


15.93 


0.14 


18 


-95.3 


8100 ± 100 


3.00 ±0.10 


2 


11.55 ± 0.63 


B218 


15.99 


0.16 


16 


-96.7 


8100 ± 200 


2.95 ±0.05 


2 


14.88 ± 0.69 


B78 


15.99 


0.15 


34 


-110.8 


8200 ± 100 


3.05 ±0.05 


2 


9.45 ± 0.25 


B153 


15.95 


0.14 


7 


-113.8 


8300 ± 200 


2.95 ±0.05 


2 


29.25 ± 1.35 


B244 


15.96 


0.14 


23 


-116.0 


8300 ± 300 


3.05 ±0.10 


1 


9.59 ± 0.41 


B331 


16.04 


0.14 


12 


-108.4 


8300 ± 200 


3.05 ±0.05 


2 


7.93 ± 0.34 


B177 


16.03 


0.15 


25 


-110.6 


8600 ± 300 


3.00 ±0.05 


2 


10.70 ± 0.37 


B424 


15.89 


0.14 


6 


-106.5 


8600 ± 200 


3.05 ±0.10 


2 


35.81 ± 3.63 


B130 


15.96 


0.15 


22 


-115.2 


9000 ± 1000 


3.05 ±0.10 


2 


5.07 ±0.24 


B267 


16.72 


0.03 


109 


-115.2 


10000 ± 100 


3.55 ±0.10 





7.22 ± 0.17 


B334 


16.58 


0.02 


8 


-109.5 


10800 ± 1000 


3.55 ±0.10 


2 


11.86 ± 1.25 


B348 


16.69 


0.01 


7 


-107.7 


11600 ± 800 


3.60 ±0.10 


2 


11.20 ±2.68 


B84 


16.56 


0.00 


49 


-109.0 


11700 ± 1000 


3.60 ±0.10 





6.12 ± 0.20 


B279 


16.56 


0.01 


144 


-104.8 


12000 ± 1000 


3.60 ±0.10 





6.25 ± 0.11 


B203 


16.68 


-0.01 


57 


-95.2 


12200 ± 1000 


3.60 ±0.10 





6.47 ±0.27 


B315 


16.80 


-0.02 


48 


-104.6 


12800 ± 800 


3.75 ±0.10 





4.26 ± 0.18 


B374 


16.79 


-0.02 


91 


-107.2 


13000 ± 1000 


3.70 ±0.10 





5.10 ± 0.16 



mean and dispersion 



-106.6 ±7.0 



